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Abstract 

The status of heat and work in nonequilibrium thermodynamics is quite confusing and non- 
unique at present with conflicting interpretations even after a long history of the first law dE(t) = 
d e Q(t) — dW c {t) in terms of exchange heat and work, and is far from settled. Moreover, the 
exchange quantities lack certain symmetry (see text). By generalizing the traditional concept to 
also include their time-dependent irreversible components d\Q(t) and d\W(t) allows us to express 
the first law in a symmetric form dE(t) = dQ(t) — dW(t) in which dQ(t) and work dW(t) appear 
on equal footing and possess the symmetry. We prove that diQ(t) = diW(t); as a consequence, 
irreversible work turns into irreversible heat. Statistical analysis in terms of microstate probabilities 
Pi(t) uniquely identifies dW(t) as isentropic and dQ{t) as isometric (see text) change in dE(t), a 
result known in equilibrium. We show that such a clear separation does not occur for d c Q{t) 
and dW c (t). Hence, our new formulation of the first law provides tremendous advantages and 
results in an extremely useful formulation of non-equilibrium thermodynamics, as we have shown 
recently [Phys. Rev. E 81, 051130 (2010); ibid 85, 041128 and 041129 (2012)]. We prove that an 
adiabatic process does not alter pi. All these results remain valid no matter how far the system is 
out of equilibrium. When the system is in internal equilibrium, dQ{t) = T(t)dS(t) in terms of the 
instantaneous temperature T(t) of the system, which is reminiscent of equilibrium, even though, 
neither d e Q(t) = T(t)d e S(t) nor d\Q{t) = T(t)d\S{t). Indeed, d[Q(t) and diS(t) have very different 
physics. We express these quantities in terms of d e pi(t) and diPi(t), and demonstrate that Pi(t) has 
a form very different from that in equilibrium. The first and second laws are no longer independent 
so that we need only one law, which is again reminiscent of equilibrium. The traditional formulas 
like the Clausius inequality § d c Q(t)/To < 0, A C W < —A [E(t — ToS(t))], etc. become equalities 
§ dQ(t) /T(t) = 0, AW = —A [E{t — T(t)S(t)], etc, a quite remarkable but unexpected result in 
view of the fact that A;5(t) > 0. We identify the uncompensated transformation N(t, r) during a 
cycle. We determine the irreversible components in two simple cases to show the usefulness of our 
approach; here, the traditional formulation is of no use. Our extension bring about a very strong 
parallel between equilibrium and non-equilibrium thermodynamics, except that one has irreversible 
entropy generation diS(t) > in the latter. 
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I. INTRODUCTION AND CONTROVERSY 



Heat and Work 



Gislason and Craig jl| recently remarked that the definition of work in nonequilibrium 
"...thermodynamics processes remains a contentious topic," a rather surprising statement, 
as the field of thermodynamics is an old discipline. However, there is some truth to their 
critique, which was motivated by an earlier paper by Bertrand [2|, who revisited the confusion 
first noted by Bauman [3] about different formulation of work [the exchange heat work 



d e W(t) = P dV(t) or the generalized work dW{t) = P(t)dV(t), see Refs. 14- 
in terms of the external (Pq) and instantaneous or internal (P(t)) pressures 
instant t, see Fig. (TJ and discussed by many others 



12 



and below] 
11| at a given 



since then with no consensus 



12j, p. 181, Vol. 



thermodynamics 



34 



Traditional formulation of nonequilibrium statistical mechanics and 



35] follows a mechanistic approach in which the system of interest 



follows its (classical or quantum) mechanical evolution in time. Being a mechanical concept, 
work is easier to identify in thermodynamics than heat; the latter is not possible to be 
identified with any mechanical force. Therefore, one usually identifies heat by first identifying 
work and then subtracting the latter from the energy change 34| in the first law. This means 
that different formulations of work will result in different heats. We refer the reader to the 

rn 

above references for an interesting history of the confusion. Gislason and Craig [29] list 
twenty-six representative textbooks including J^f], 0] where the pressure-volume work and 
heat are defined so differently that they are not equivalent in the presence of irreversibility, 



and there appears to 
debate in the field 
variables 



3e no consensus about their right formulation so far; see the recent 



17 



21 



7. 8. 



38 



26(| . Unfortunately, none of the latter sources consider internal 



and references therein] that are needed to describe irreversible processes 
but do not appear explicitly in the Hamiltonian. They also do not consider thermodynamic 
forces js S] that are non-zero when the system is away from equilibrium. 

The attempt to define heat by first defining work leaves the concept of heat devoid of 
clear physical significance, especially since it depends on what we mean by work. In many 
cases, work is identified by considering the work performed by bodies external to the system 
(the medium; see Fig. [T|), which may have nothing to do with the work done by the system, 
especially when we consider irreversible processes. This has created a lot of confusion in the 
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literature when dealing with irreversible processes; see Refs 



sion; see also 



34 



26 1 for a recent discus- 
10| . The traditional formulation of statistical mechanics and thermodynamics 



35] follows a mechanistic approach in which the system of interest, to be precise its mi- 



crostate, follows its (classical or quantum) mechanical evolution in time. The only difference 
between mechanics and thermodynamics is that the evolution in thermodynamics is always 



stochastic, which makes t 
example a recent review 



re evolution irreversible in accordance with the second law; see for 



39] . Therefore, it is crucial to consider a statistical foundation of 



" , " i " 1,,m ™*^*^™^»^"» " 



have also initiated recently |7H9l. l39l. |40||. The concept of work follows in a trivial manner by 
taking (a stochastic average of) its mechanical analog (work done by generalized forces in 
the Hamiltonian formulation) [34j. and was shown to be given by dWy(t) = P(t)dV(t) for 



the pressure work; see Ref. 



7] . These forces can be controlled by an observer and determine 



the observables in the system; see for example sj]. It was shown in Ref. 7j,|8| that the first 
law can also be written in terms of the generalized heat dQ(t) added to and the generalized 
pressure work dW{t) done by the system: 



dE(t) = dQ(t)-dW(t), 



(1) 



where 



diQ(t) = dQ{t) - d e Q{t), diW{t) = dW{t) - d c W(t), (2) 

denote irriversible heat and work, respectively, generated within the system {4]. However, the 
interpretation of the above new formulation, its statistical basis, and its realtionship with the 
traditional formulation dE(t) = d c Q(t) — d c W(t) in terms of exchange hea t dp Q(t) = T d e S(t) 
and work d e W(t) = P dV(t) was not explored in earlier publications [7- 10] . 



For an isolated system, the mechanical observables, collectively denoted by the set X, 



are additive integrals of motion 



34 



41] such as the energy E, the number of particles Nj 



of different species, and linear and angular momenta P and M of the system, respectively, 
whose values depend on the preparation of the system. When the motion is confined to a 
finite region of space, the volume V of this region also characterizes the system as a con- 
straint. The constraints are also treated as observables. The internal variables unfortunately 
cannot be controlled by an observer. Despite this, the concept of work needs to incorporate 
the additional work required in changing the internal variables. This gives rise to some 
complication in identifying the form of work associated with them. The collection of ob- 



servables and internal variables will be called state variables and will be denoted by the 
set Z. The macrostate can be characterized in a larger state space of Z, even though the 
observed macrostate is characterized by X. The internal variables do not remain constant 
for an isolated system as the system relaxes. In this sense, we must treat the isolated system 
as interacting with an appropriate "fictitious" medium to allow for the variation of internal 
variables, a fact not emphasized in the literature to the best of our knowledge; see Ref. Q]. 
As both energy and work are mechanical concepts, it is trivial to identify their statistical 
averages over microstates whether the system is in or out of equilibrium. 



B. Controversy 

As the concept of heat has no analog in mechanics, it is usually identified indirectly as that 
energy change which is not work. This kind of approach makes heat and work not unique: 
a change in one affects the other; only their difference has a well-defined meaning. This has 
created a lot of confusion in the field. Our goal in this work is to establish that there exists 
a consistent procedure in which the partition is unique so that work and heat are no longer 
arbitrary in irreversible processes. They are given by generaized work and heat despite the 
controversy. The first law of thermodynamics, see for example, Refs. @, Q]> * s a ver y general 
statement and is supposed to be valid for all processes, and not just equilibrium processes. 
The traditional formulation of the first law in terms of the exchange heat and work is not 
only oblivious to the internal variables and thermodynamic forces, whose presence must 
control the approach to equilibrium, but most significantly, is also oblivious to the violation 
of the second law; see below. This is why we need to use both the first and the second laws 
in nonequilibrium thermodynamics. How can a formulation of a fundamental law of physics 
allow for the violation of another fundamental law of physics? There is also an asymmetry 
between heat and work that will be elaborated below in that the two terms are not on an 
equal footing. There is obviously no problem for reversible processes. 

We mostly consider mechanical work including dissipation, but the arguments are valid 



for all kinds of work; see Sec. IVIIIt however. Zemansky 37], p. 73] defines heat as energy 
exchange "...by virtue of a temperature difference only." Unfortunately, this rules out any 
isothermal heat exchange and cannot be considered general. Kirkwood and Oppenheim de- 
fine heat as energy exchange resulting in "...the temperature increment..." (which rules out 
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phase changes requiring latent 



36 



reat) and later note that the work may be converted to heat 
pp. 16,17] as was first observed by Count Rumford in 1798 
. These are two of the examples where heat is defined directly 



due to frictional dissipation 
during the boring of cannon 
without first identifying work; they have limitations and cannot be considered general, es- 
pecially when irreversibility is present. The situation with work is just as confusing as we 
noted above with its two different formulations: d e W(t) = PodV(t) or dW(t) = P(t)dV(t). 

Dissipation always gives rise to positive entropy generation due to irreversibility and also 
raises the temperature such as due to friction or the Joule heat in resistors. Therefore, 
it is natural to account for such viscous dissipation in work and heat when dealing with 
nonequilibrium systems, as they are integral to the system and dictate its relaxation. The 



fact that literature is not very c 
this work; see however 



31 



to obtain the results derived here. Recent 
in stochastic trajectory thermodynamics 



ear on how to incorporate viscous dissipation has motivated 
33(, but the authors do not take the discussion far enough 
y dissipative forces are explicitly considered 



15 



17- 



our approac 



li 



271 ] . However, the approach differs from 



in important ways and has also given rise to controversy, and remains 



contentious 17], l2lM26j . In addition, the approach is not general as it is limited to isothermal 
variations and to cases where Langevin dynamics is applicable; see, however, Hoover and 
Hoover |43j where an example of a time-reversible deterministic Hamiltonian system is given. 



C. New Results 

We interpret heat and work used in the recent reformulation in Eq. (TjQ) that was proposed 
in Refs. 7H9| and follow the consequences with an aim to develop their statistical definition. 
The formulation makes the first law identical to the second law (in the guise of the Gibbs 
fundamental relation) In general, d c Q(t) = T d e S(t), but dS(t) and dQ{t) are not related to 



each other, except in internal equilibrium [7[. Moreover, it would be incorrect to conclude 
from dS(t) = d e Q(t)/T + d[S(t) that d[Q(t) = T^diS^); see later. The aim of any theory of 
nonequilibrium thermodynamics is to determine the entropy change dS(t). Therefore, the 
determination of diS(t) becomes the focus in any investigation of a body in the traditional 
formulation. The second law is reflected in the inequality diS(t) > for any irreversible 
process. The statistical analysis provides an elegant formulation of nonequilibrium thermo- 
dynamics in which there is a unique and natural distinction between nonequilibrium work 
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and heat in that the generalized work represents isentropic change in the (internal) energy 
and the generalized heat represents the isometric (constant extensive state variables exclud- 
ing the energy) change in the energy. The latter change results purely from the entropy 
change. The unique partition of energy remains valid no matter how far the system is out 
of equilibrium. The assumption of internal equilibrium allows us to express the first law 
in terms of the instantaneous (or internal) fields. This brings about a very close parallel 
between nonequilibrium and equilibrium processes such as dS(t) = dQ(t)/T(t) and the ex- 
istence of a theorem for irreversible processes identical in spirit to the adiabtic theorem 34] 
for equilibrium processes. The approach proposed recently in Refs. UB Q deals directly 
with dS(t) without a need to use d e S(t) and diS(t) separately, although they can also be 
evaluated in our approach. Therefore, our approach should be quite useful as the entropy is 
a state variable; see below, however, for other advantages. A particular symmetry is explic- 
itly seen in our formulation in that both dQ and dW not only do not change with the nature 
of the process but also exhibit an identical formulation in terms of entropy and volume, 
respectively. In other words, they are found to be on an equal footing. As a consequence, 
the first law becomes identical to the second law (as the Gibbs fundamental relation) so we 
only deal with a single law. The Clausius inequality (4^] turns into an equality in all cases 
as dQ(t)/T(t) becomes a state variable, the work is expresssed as an equality, as if we are 
dealing with equilibrium processes, a quite remarkable result in its own right, even though 
there is irreversible entropy generation. It has been recently suggested [32] that the use of 
internal fields is not always consistent with the second law. We find no such problem in our 
approach. 

The new formulation of the first law contains not only the observables and the instanta- 
neous fields of the system but also explicitly contains internal variables so no information 
about the system is lost. This makes the new formulation quite beneficial. The instan- 
taneous temperature, pressure, etc. of the system need not be identical to those of the 
medium when the system is out of equilibrium with the medium as has become clear re- 
cently Q UR Thus, o Ur approach a. 8 o overcomes the object.on raM by Cohen and 



Mauzerall [21J against the stochastic trajectory thermodynamics 15|, [l7|, |27J to which our 
approach is easily extended. We will establish that there are several benefits in accounting 
for viscous dissipation within the system in the first law, but care must be exercised. Our 



approach, which assumes the existence of entropy even for nonequilibrium states 



39] as 
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FIG. 1: Schematic representation of S, S and So- We assume that S and S are homogeneous and 
in internal equilibrium, but not in equilibrium with each other. The internal fields T(t),P(t), ■ ■ ■ 
fof £ and Tq,Pq, - ■ ■ of £ are not the same unless they are in equilibrium with each other. There 
will be viscous dissipation in £ when not in equilibrium with £. 



was first proposed by Clausius [44], is otherwise very general as we do not restrict ourselves 
to any particular dynamics and to only isothermal variations. It allows us to incorporate 
viscous dissipation, due to the presence of thermodynamic forces and internal variables, in 
the discussion explicitly within this general framework. We consider our system E (see Fig. 
[Q surrounded by a very large medium E so large that its fields such as its temperature 
T , pressure Pq, etc. are not affected by the system. They form an isolated system S . We 
consider all systems to be stationary; their relative motion will be considered only in Sec. 
IVHI The surface separating E and E, which may represent a piston, will be treated as having 
no interesting thermodynamics of its own although it may participate in irreversibility due 
to field (such as temperature, pressure, etc.) differences across it. Because of the enormous 
size of the medium with respect to that of the system, all irreversible components in E 
appear within E (d\Q{t) = diQo(t),diW(t) = diWoit)) although this is not always so in the 
literature 46]. The situation of finite surroundings has been considered by Bizarro [31] and 
by us in Ref. jsj].) In the following, all extensive quantities pertaining to E and So carry 
an annotation tilde or a suffix 0, respectively, and those pertaining to E carry no suffix. We 
will use body to refer to any one of the above three systems and use symbols without any 
suffix to denote its quantities. 

Our discussion is easily extended to the case when E is comparable to E in size, as is easily 
seen in Sees. IVIII and IIX | and Ref. [8j. In this case, we will refer to E as the surroundings, 



S 



whose fields will also change with time. However, to keep the discussion simple, we will 
treat S as an extensively large medium for the most part. 

The layout of the paper is as follows. We briefly review the traditional formulation of 
the first law, its lack of symmetry, confusion and its limitations in the next section. In 
Sec. 1111} we follow our newly developed nonequilibrium thermodynamics and unravel the 
significance of heat and work in that approach. Several important results are derived there 



in the forms of theorems. Some of the results were announced earlier in Ref. but have 
now been expanded. In Sec. IIV} we argue that only the generalized work is consistent with 
the second law. The statistical definition of heat and work is taken up in Sec. |V] which 
forms the core of the present work, where we show that this definition is identical with heat 
and work discussed in Sec. IIIII and used in the earlier work |z-9]. This section also contains 
many important results including the adiabatic theorem for irreversible processes. A general 
expression for microstate probabilities is derived here, which clearly shows nonequilibrium 
effects in their formulation. We partition microstate probability changes into external and 
internal parts d e pi(t), d\Pi(t) to obtain the statistical formulations for d e Q(t), d\Q{t) etc. We 
use the example of an ideal quantum gas to show how d e pi(t) , d\pi(t) can be computed. In 
the following section, we discuss the Clausius inequalities (there are two different ones) and 
the Clausius equality and the work equality. In Sec. IVIfl we consider two applications of 
our formulation, where the traditional formulation cannot be applied. In Sec. IVHH we 
extend our discussion to include an additional observable. The results derived here are used 
in Sec. HXI to study a closed system which is allowed to exchange some kind of "work" with 
a thermally isolated external object. This is a classic prototype model studied extensively; 



see for example, Ref. 34j . The last section contains a brief summary of results and a list of 



benefits of our approach. 

II. TRADITIONAL FORMULATION OF THE FIRST LAW 
A. Traditional Formulation 

To truly appreciate our contribution, it is useful to consider how the first law is tradition- 
ally expressed. We will only consider a single internal variable £ for simplicity. Similarly, we 
will consider only E and V for simplicity as observables with the number of particles N (only 
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a single species) held fixed, unless noted otherwise. Other variables are easy to include in the 
approach as we do in Sec. IVIIII Traditionally, dQ{t) represents the amount of heat exchange 
d e Q(t), so that —dW{t) is identified with the work exchange —d e W(t) = dW(t) = P dV(t) 
by the medium to the system, giving d e W(t) = PodV(t). This is true even if the system 
cannot be assigned any pressure or if its instantaneous pressure P(t) is different from P Q . As 
the net heat exchange d c Q(t) + d c Q{t) = 0, we immediately verify d c Q{t) = T d c S(t). There 
is no such general relation relations for other heats: dQ(t) ^ TodS(t) and d\Q(i) ^ T^d^Sit). 
The traditional formulation of the first law for a general process reads 

dE(t) = d e Q(t) - d c W(t) ee T d c S(t) - P dV(t) (3) 

expressed in terms of either exchange quantities or external fields of E. Only when the pro- 
cess is reversible that we have dE(t) = TodS(t) — PodV(t). The external fields are conjugates 
to the observables in X (E, V), with the medium affinity Aq conjugate to £ vanishing. Thus, 
the above formulation of the first law is, as said earlier, oblivious to the internal variables 
and will not be considered when using this formulation. The following inequality for a cycle, 
commonly known as the Clausius inequality, follows from dS > d c Q(t)/T , 

§d c S(t) ee §d c Q{t)/T < 0, (4) 



B. Confusion about Work and Heat 



The situation regarding dW(t) is not always clear. Kondepudi and Prigogine use dW(t) = 
PdV(t), where P " ...is the pressure at the moving surface," but they do not mention whether 
the form is applicable to all processes. Landau and Lifshitz are explicit and state that 
dW(t) = P(t)dV{t) for reversible and irreversible processes 34j, p .45j. They require for this 



the existence of mechanical equilibrium (and so do Refs. [36|, l37|) within E so that at each 
instant during the process P(t) must be uniform throughout the body; its equality with P 
is not required. However, they do not discuss dQ(t) when they consider E out of equilibrium 



with E 



3J, Sect. 20]. If we use dW(t) = P(t)dV(t) for the work done by E, then this will 



alter the heat dQ(t) added to E in Eq. ([T]). This follows immediately from the fact, not 
appreciated in the literature to the best of our knowledge, that dE(t) must be invariant to 
the choice of internal or external fields. Also, to the best of our knowledge, the issue of the 
actual forms of dQ{t) and what is the correct form of dW(t) for nonequilibrium processes 
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has not been settled in the literature. Indeed, Kestin [12|, Sect. 5.12] clearly states that 
distinguishing heat and work in nonequilibrium states is not unambiguous. We will later 
argue otherwise in this work. 



C. Lack of Symmetry and Disconnection with the Second Law 

The exchange heat and work d c Q(t) and d e W(t) are very different for irreversible pro- 
cesses, since d e Q{t) = T d c S(t) is in terms of d e S(t), which is not the change in a state 
variable, while d e W(t) = P dV(t) depends on the change in a state variable. Thus, there 
is an asymmetry between the two in the second equation in Eq. ([3]) in that they are not 
on an equal footing. However, the main disadvantage of the formulation is that it is always 
valid, even for a process that violates the second law by having the exchange heat flow from 
a colder to a hotter object. This is why we need both laws in the traditional formulation of 
non-equilibrium thermodynamics. We believe that an elegant formulation of a fundamental 
law like the first law should not only satisfy the other fundamental law, the second law but 
also exhibit as much symmetry as possible. 

D. Limitations 

As the Gibbs fundamental relation, see Eq. ([7]), for the system explicitly contains the 
internal variable, the first law in its traditional formulation and the Gibbs fundamental 
relation, which codifies the second law, have different contents for irreversible processes; 
in particular, the latter contains more information than the former. Only for reversible 
processes, for which internal variables are no longer independent of the observables, the first 
law can be used to determine the change AS(t) = [dE(t) + PodV(t)]/To in the entropy of a 
body. This is not true when we deal with irreversible processes. This limits the usefulness 
of the first law. We now list some of the important limitations of the traditional formulation 
below. 

(1) The law is oblivious to the violation of the second law. 

(2) There is the above mentioned asymmetry between heat and work. 
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(3) As both d c Q(t) and d e W(t) are determined by the medium, their knowledge does not 
provide us with any direct information about the system or its entropy change AS(t) 
in an irreversible process. This is easily seen by considering an isolated body. As 
d c Q(t) = d c W(t) = 0, the first law has no useful content. Moreover, it cannot be used 
to determine dS(t) = rfiS'(t) as the body changes from some state A' to a nearby state 
A unless both states are equilibrium states. In the latter case, one can compute dS(t) 
by consider some equilibrium path connecting the two states. This approach will not 
work if one of the two states is or both are out of equilibrium. 

(4) Work and heat cannot always be unambiguously distinguished, a point already made 
very strongly by Kestin [121 , Sect. 5.12] and which is at the heart of the dispute 
discussed above. 



(5) The heat and work do not always flow through the boundary 12|, footnote on p. 176]. 



(6) As d e S(t) can be determined from d e Q(t), the aim of any nonequilibrium thermody- 
namic investigation using the traditional formulation is to determine the irreversible 
entropy change diS(t). For this, one needs to invoke the Gibbs fundamental relation 
in addition to the traditional first law; see for example de Groot and Mazur j^J. 



(7) As internal variables £(t) |5j,|7H9 



38 



, |45|, |47J , which are very common in nonequilibrium 



systems such as glasses or in chemical reactions, play an important role in nonequi- 
librium thermodynamics, their behavior will strongly affect the dissipation within the 



system. Unfortunately, these variables do not couple to the medium 38]; hence, they 
do not appear in d e Q(t) and d e W(t), although they control the thermodynamic relax- 
ation and, in particular, the Gibbs fundamental relation for the system. 

III. GENERAL CONSIDERATION 

As a generalization of the equilibrium concept, we use the instantaneous values of the 
state variables containing the observables (X(t) :E(t),V(t)) and internal variables (£(£)) 



t| to identify the state of the body. A body can be in the same state at different 
times. This is important so that a system can go through a cyclic process in which the 
system comes back to the same initial state at a later instant. However, the entropy of 
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a body at some instant, besides being a function of the state variables, may also have 
an explicit dependence on time: S(t, E(t), V(t), £(£)). Thus, in general, the entropy will 
not be a state function. The first and second laws are not useful for any computation 
unless we can ascribe temperatures, pressures, etc. to S. This requires S and £ to be in 
internal equilibrium {7, ^ [34 1 when their instantaneous entropies become state functions 
S(t) = S [E(t),V(t),£(t)] ,S(t) = S E(t),V(t),£(t) of (time-dependent) state variables. 
Let W{t) = W[E(t),V(t),£(t)] and W(t) = W \E(t),V(t),£(t)\ denote the number of 
microstates consistent with the state variables for each of them. Then, as discussed in Refs. 
[7, 8], we have 

S(t) = hxW(t),S(t) = lnW{t) (5a) 

in other words, the microstates in W(t) or in W(t) are equally probable. 

The temperatures, pressures and affinities (we introduce (3(t) = 1/T(t) and (3q ee l/T ) 
are given by appropriate standard derivatives of the entropies: 

P{t) = dS(t)/dE(t), /3(t)P(t) = dSit)/8Vit), Pit) Ait) = (dS(t)/d£(t); (6a) 
0o = dSit)/dEit), (3 P = 8Sit)/dVit), dSit)/dVit) = 0. (6b) 

The Gibbs fundamental relations are given by 



dE(t) = T(t)dS(t) - P(t)dV(t) - A(t)d£(t), dE(t) = T dSit) - P dV{t); (7) 

The validity of Eq. ([7| requires S and S to be independently homogeneous 3, |8j and in 
internal equilibrium. We now prove the following trivial but important theorem. 

Theorem 1 Irreversible work and irreversible heat have identical values: 

diQ(t) = diW(fi), (8) 



Proof. We have dEo = dVo = Aq = for So- The application of the first law for £0 using 
generalized heat and work, see Eq. fl2]), gives 



dE = dQ it) - dWoit) = diQo(t) - diW (t) = 0; 



(9) 



there is no exchange heat and no exchange work for So- As the irreversibility is only 
associated with the system E, we have d\Q{t) = d^Q^it) and d\W{t) = diWo(t). The desired 
equality in Eq. (jSJ) now follows from Eq. (jUJ). ■ 
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The statistical demonstration of the identity in Eq. (JSJ) is given in Theorem [6j It is 
now easy to show that heat and work can also be used in the first law for the system 
(dE(t) = dQ(t) - dW(t)), since dQ{t) - dW{t) = d e Q(t) - d e W(t). For the medium, it also 
holds (dE(t) = dQ(t) — dW(t)) as d\Q(i) = d\W(t) = 0. Thus, we can express the first law 
for any body by also using generalized heat and work; see Eq. (JTJ). 

It was established in Ref. Q that 

dQ(t) = T{t)dS{t). (10) 

It is a generalization of d e Q(t) = T Q d e S(t) to dQ(t); the latter denotes the heat added to the 
system either through exchange with its exterior (d e Q(t)) or by dissipative internal forces 
within (diQ(t)). Similarly, 

dW(t) = d c W(t) + ckW(t) = P(t)dV(t) + A(t)d£(t) (11) 

is the generalization of work done by the system: it includes the work done on its exterior 
(d e W(t)) and the (internal) work done by dissipative internal forces (diW(t)). 

To appreciate the importance of the new definition of heat and work, let us for the 
moment assume that there is no internal variable and that P(t) > Pq. The irreversible work 
diWo(t) done in E by the pressure difference AP(t) = P(t) — Pq > is 

dW (t) = diW (t) = AP(t)dV(t) > 0, 

since dV(t) > 0, and appears as the irreversible work within the system and results in 
raising the kinetic energy dK$ of the center-of-mass of the surface separating S and S and 
overcoming work dWf T (t) done by all sorts of viscous or frictional drag. Thus, 

diW (t) = dK s + dW^t). (12) 

Because of the stochasticity associated with any statistical system, both energies on the right 
side dissipate among the particles so as to increase the entropy and appear in the form of heat 



(diQ (t) = diQ(t) > 0) within the isolated system |39|. Thus, when there are irreversible 
processes going on, it is natural to generalize heat from d e Q(t) in Eq. ([3]) to include d^Qit) = 
d\Qo{t) and identify dQ(t) as the heat added to the system. Similarly, we need to generalize 
work from d e W(t) = PodV(t) to dW(t) = P(t)dV(t) and identify it as work done by the 
system. In the presence of the internal variable, there is an additional contribution (A(t) — 
A )dt;(t) = A(t)dt;(t) > to dW(t). This does not change the conclusions above. 
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We finally conclude that 

dE(t) = d c Q(t) - d c W(t) = dQ(t) - dW{t) (13) 

which demonstrates that both formulations are valid for the first law. However, the most 
important result is given by Eq. (flOl) . We also see that d\Q{t) ^ T diS(t), d{Q(t) ^ 
T(t)diS(t), even though d e Q(t) = T d c S(t); see Eq. (fTBj) . 

The equality in Eq. (jSJ) can also be obtained by the use of the Gibbs fundamental relation. 
We follow the approach initiated in Ref. {sj], and rewrite dE(t) by explicitly exhibiting the 
thermodynamic forces as follows: 

dE(t) = T d c S(t) - P dV(t) + T diS(t) + [T(t) - T }dS(t) + [P - P(t)}dV(t) - A(t)d£(t) 
to conclude that 

T diS(t) + [T(t) - T ]dS(t) + [P - P{t))dV(t) - A(t)d£(t) = 0. (14) 

For this to be valid, each of the last three terms must be non-positive: 

[T - T(t)]dS(t) > 0, [P(t) - P }dV(t) > 0, A(t)d£(t) > 0, (15) 

to ensure that d[S(t) > 0. The factors T Q —T(t), P(t) — Pq and A(t) in front of the extensive 
state variables are the corresponding thermodynamic forces that act to bring the system to 
equilibrium. In the process, each force generates its own irreversible entropy generation Q]. 
The equalities occur when thermodynamic forces vanish. 

It is useful to acknowledge at this point that there are no thermodynamic forces in the 
medium. To see this, we consider Eq. ([7j) for dE(t), in which both terms contain the constant 
fields Tq and Pq of the medium, clearly showing that the thermodynamic forces are zero. 
This means that there cannot be any irreversible entropy generation within the medium; 
they only appear within the system, as we have said earlier. 

We now recognize that 

diQ(t) = Tod.Sit) + [T{t) - T )dS{t), (16a) 
= T(t)<kS(t) + [T(t) - T ]d c S{t) (16b) 

and 

ckW(t) = [P(t) - P ]dV(t) + A(t)d£{t). (17) 
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Their equality merely reflects the fact that in the partition dE(t) = d e E(t)+diE(t), d\E(t) = 
0. We also note that while each term in d\W(t) is non-negative, this is not so for d\Q(t) in 
which the first term is non-negative, but the second term in Eq. ( II 6p is non-positive. This 
not only means that the physics of diQ{t) and d[S(t) is very different but also that 

diQ(t) < T diS(t) , dQ(t) < T dS(t); (18) 

the equalities occur only for isothermal (T = To) or adiabatic (dS = 0) processes. Even 
though dQ(t) = T(t)dS(t) for a system not in equilibrium with the medium, we have 
d e Q(t) T(t)d e S(t), d\Q(t) ^ T(t)d\S(t). It should become evident by now that it would be 
incorrect to use dS{t) = d c Q(t)/T + diS(t) to conclude d- x Q(t) = T Q diS(t). 
In the general case, the first law can be written as 

dE(t) = S^-dZ(t), 
V ; dZ(t) K 1 

where Z(t) contains S(t) and the set Z#(t) consisting of all state variables except E(t). We 
clearly see that each term in the scalar product has the same mathematical form, ensuring 
that all terms are on an equal footing. It is this symmetry that was absent in the traditional 
formulation, but is present in the new formulation of the first law. As this is also the general 
form of the Gibbs fundamental relation, the two laws have reduced to a single law, as we 
have claimed. Therefore, out formulation of the first law will always remain consistent with 
the second law. 



IV. GENERALIZED WORK AND THE SECOND LAW 

Let us follow the consequences of this particular generalization a bit further by 
again restricting to no internal variable for simplicity, and prove that only dW(t) = 
P(t)dV(t),dW(t) = P dV(t) using the internal pressures of the bodies is consistent with 
the second law, and not dW(t) = PodV(t),dW(t) = P(t)dV(t) which use the pressures 
external to the bodies. These choices for work are symmetric as opposed to the traditional 
formulation in Sec. [Til in which there is no symmetry between d e W(t) and d c W(t). 

We take P(t) > P and consider the choice dW(t) = P(t)dV(t) etc. We obtain dW(t) + 
dW(t) = P(t)dV(t) + P dV(t) valid for any arbitrary dV(t) = -dV(t) so that 

diW{t) = [P(t) - P ] dV(t) > 0, (19) 
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which is consistent with the second law and proves the above assertion, once we recognize 
that dV(t) > 0. The second choice will result in the violation the second law, since dW\(t) = 
(P — P{t))dV(t) < 0, a physical impossibility. Thus, we must write the first law for the sys- 
tem and the medium, respectively, as dE(t) = dQ(t) — P(t)dV(t), dE(t) = dQ(t) — P dV(t). 
The above discussion is easily extended to include internal variables without affecting the 
above conclusion. 

The above demonstration establishes that the work done by a body is given by Eq. 
(II ip in all cases contrary to the traditional formulation, see Eq, fl3]), in which it is given 
by P dV(t). The generalized formulation brings out the another symmetry: under the 
interchange system<^=^medium, work and heat for any body always uses its own internal 
fields. This symmetry is absent in the traditional formulation. The new symmetry will 
prove very useful when the medium is not extensively large compared to the system or when 
we need to consider mixing of gases, free expansion, etc. where there is no clear separation 
between different parts of an isolated system into a system and a medium. 

V. STATISTICAL DEFINITION OF WORK AND HEAT 
A. System 

1. System not in Internal Equilibrium 

Before proceeding further, let us see how the generalized heat and work could be under- 
stood from a statistical point of view. We consider two possible neighboring nonequilibrium 
states A and A' at different times t and if < t, respectively, so that the differences in their 
state variables dE(t) = E(t) - E'(t'), dV(t) = V(t) - V'(t'),d£(t) = £(t) - f'(f) and the 
difference dS(t) = S(t) — S'(t') in their entropies are infinitesimal. We use the index % to 
label the microstates of the system and let Pi(t) , p'i(t') their their probabilities in A and A', 
respectively. These probabilities are functions of the state variables (including the number of 
particles, but that remains constant) and may also have an explicit time dependence. Thus, 
the discussion here does not require the system to be in internal equilibrium. Obviously {48 1 

5>(*) = Edit) = 1. (20) 
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We will see below, see Theorem [6j that the probability conservation is behind the statistical 
demonstration of the identity in Eq. (jHJ) . The energy Ei of the ith microstate, on the other 
hand, depends on V(t) and [in general, Ei will depend on the set Z E (t)}, but will have no 
explicit t-dependence. The entropy S and the energy E are given by the following averages 

S(t) ee EMtMt), E(t) = EMt)m), (21) 

where 

^EE-ln^(t) (22) 

is the uncertainty of Shanon or the negative of the index of probability of Gibbs (9), [soj] and 
Ei(t) is the energy of the ith microstate. (We will avoid the use of microstate "entropy," to 
refer to r/j, which has become common in the literature these days.) The entropy expression is 



due to Gibbs 35]. We have exhibited a time-dependence in Ei(t) to reflect the fact that this 
energy can change as V(t) and £(£) change during the transition A' — >A. The microstate 
probability will also change in time. In particular, a microstate may disappear or a new 
microstate may emerge in time. This is most easily seen by recognizing, see Eq. ([5]), that 
the entropy is determined by W(t). As entropy changes, W(t) must change so that either 
some previous microstates disappear or some new microstate emerge. All this will become 
clear below in Sec. IV Gl where we discuss a simple example of an ideal gas. We now prove 

Theorem 2 E(t) is a function of V(t),£(t) and S(t), even though Ei[V(t), £(t)] are func- 
tions ofV(t) and £(t) only. 

Proof. We consider the differential 

dE(t) = ZiEi{t)d Pl {t) + Y.iPi(t)dEi(t). 

As Pi(t) are unchanged in the first sum, this sum is evaluated at constant entropy. Thus, 
this contribution is isentropic which we denote by dE\ s . The microstate energies E{ are 
unchanged in the second sum so this contribution refers to an isometric process at fixed 
V{t) and ^(t) and we denote the contribution by dE\ v ^. Thus, 

dE = dE\ v ^ + dE\ s . (23) 

This proves that E(t) is a function of S(t),V(t) and ^(t). ■ 
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In general, E(t) is a function of S(t) and the set Z^(t). We introduce a special process, to 
be called an isometric process, which is a process at constant Z#(t) and is a generalization 
of an isochoric process. In this process, the work done by each mechanical variables in Z#(i) 
remains zero. We now prove the following theorem that establishes the physical significance 
of the two contributions. 

Theorem 3 The isentropic contribution represents the generalized work dW(t) and the iso- 
metric contribution represents the generalized heat dQ(t). 



Proof. We follow Landau and Lifshitz 34| and rewrite the second term in Eq. (|23|) as 



dE\ s = EMt)-^dV(t) + j2. Pi (t)-^dt(t) = -J2 iPi (W(t)dv(t) - EMt)Mt)dm 

where we have introduced Pi(t) = —dEi(t)/dV(t) as the pressure produced by the ith 



microstate of the system |49|, p. 67] on its boundary. This pressure corresponds to a force 
pointing away from the system. Similarly, Ai(t) = —dEi(t)/dC,(t) as the affinity of the ith 
microstate. We assume that the changes dV(t) and d£(t) are the same for all microstates so 
that they can be taken out of the summations. In terms of 

P(t) = HiPm p i{t), A(t) = ZMVMt), (24) 

respectively, which define the instantaneous average pressure P(t) and affinity A(t) of the 
system, respectively, we can relate dE\ s with dW(t) and nit with d e W(t): 

dW(t) = - dE\ s = P(t)dV(t) + A(t)d£(t). (25) 

This identification then also proves that the heat in the first law must be properly identified 
with dQ and not with d e Q. Accordingly, 

dQ = dE\ VA ee J2Ad Pi , (26) 

i.e., dQ for irreversible processes is the isometric change in the energy, m 

We should point out that by assuming dV(t) and d£(t) above to be the same for all 
microstates, the statistical nature of dE\ s is reflected in the statistical nature of P(t) and 
A(t), the internal fields of the system. Thus, the internal fields are fluctuating quantities 
from microstate to microstate when dV(t) and d£(t) are not treated statistically. 
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In general, dE\ s will be a sum of various works dWz{t) = —(dE/dZ)s.z' E dZ, with 7i' E (t) 
consisting of all state variables in Z#(i) except Z(t) used in the derivative: 

dE\ s = -dW(t) = - dW z (t). (27) 

Z£Z E (t) 

This is again consistent with the previously mentioned symmetry in the new formulation. 

The above discussion proves that the definition of heat and work does not require the 
establishment of the internal equilibrium within the system. It is useful to compare the above 
approach with the traditional formulation of the first law in terms of d c Q(t) and d e W(t): 
both formulations are valid in all cases. It should be mentioned that the above identification 
is well known in equilibrium statistical mechanics, but its extension to irreversible processes 
and our interpreation is, to the best of our knowledge, is novel. 

We now prove a trivial but conceptually an important theorem. 

Theorem 4 Heat and work in the traditional formulation of the first law do not have a 
clear division as in Eq. (H3|). 

Proof. To prove the theorem, we focus on d e W(t) = P dV(t) and recognize that P is a 
constant. Hence, the statistical nature of d e W(t) must be contained in dV(t). As V(t) = 
J2iPiVi, we have 

d e W(t) = PoiEiPidVi + EiVidPi). 

We observe that d e W(t) contains not only an isentropic contribution, the first sum on the 
right, but also contains a contribution containing dp^. Thus, the clear separation between 
heat and work, as appears in Theorem [2], is not present in the traditional formulation of the 
first law. ■ 

2. System in Internal Equilibrium 

It should be clear from the existence of non-zero thermodynamic forces for irreversibility, 
that P(t) 7^ Pq = ( dE/dV) except when a mechanical equilibrium exists. While the 
instantaneous average pressure is defined under all circumstances, it can only be identified 
with the thermodynamic definition of the instantaneous pressure 



P(t) = - (dE/dV) Sji 
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(28) 



when the system is in internal equilibrium. Similarly, the instantaneous average affinity 
A(t) of the state has nothing to do with the affinity of the medium (Aq = 0), and can only 
be identified with its thermodynamic definition A(t) = —{dE/d^)s.v when the system is in 
internal equilibrium. To proceed further, we need to impose the condition of internal equi- 
librium, so that Pi has no explicit time-dependence. In this case, not only the instantaneous 
pressure satisfies Eq. (I28I) . but we also have, following Eq. (|26|) . 



which is a statistical proof of the thermodynamic identity in Eq. ( TTUj) relating dQ(t) and 
dS(t). We also note that the ratio dQ(t)/dS(t) is related to a field variable of a (macro)state, 
the instantaneous temperature of the system, while in the conventional approach, the ratio 
d e Q(t)/d e S(t) = T does not give a field variable of the state. 

It is clear from the above discussion that it is heat and not work that causes Pi(t), and 
therefore the entropy, to change without changing Ei. This is the essence of the common 
wisdom that heat is random motion. But we now have a mathematical definition: heat is 
the isometric part of dE(t) that is directly related to the change in the entropy through 
changes in Pi(t). Work is that part of the energy change caused by isentropic variations in 
the "mechanical" state variables (/!;). Thus, work causes Ei to change without changing 
Pi(t). This is true no matter how far the system is from equilibrium. Thus, our formulation 
of the first law and the identification of the two terms is the most general one, and applicable 
in all cases. The relationship between heat and entropy becomes simple only when the system 
is also in internal equilibrium in which case T(t) has a thermodynamic significance; see Eq. 
I[6a]) and we have the thermodynamic identity in Eq. I[W\) relating dQ(t) and dS(t). 

B. Microstate probabilities 

In internal equilibrium, the entropy must be at its maximum at fixed E(t) = 
J2i EiPi, V(t) = J2i V iPi and = J2i iiPii and is obtained by varying p t without changing 
the microstates, i.e. Ei, Vi and This variation has nothing to do with dpi in a physical 
process. Using the Lagrange multiplier technique, it is easy to show that the condition for 
this in terms of four Lagrange multipliers whose definitions are obvious is 



dQ(t)/dS(t) = (dE/dS) 



T(t), 



(29) 



rji = Ai + X 2 Ei + X 3 Vi + A 4 & 
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(30) 
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from which follows S = X\ + X2E + X 3 V + It is now easy to identify A2 = 0, A3 = f3P 
and A4 = fiA so we finally have 



Pi (t) = exp[(3(t)(G(t) -Et- P(t)Vi - A(t)£i)), (31) 
where Ai = f3{t)G{t) with G(t) defined by 

exp(-/3(t)G(t)) = J2eM-P(t)(Ei + P(W + Am)}- 

i 

The quantity G(t) would represent the free energy of the system, had it been in a medium 
E(T, P, A). However, as the system is in a medium E(T , Po, A = 0), G(t) does not represent 
the free energy in this case; the correct free energy of the system is the Gibbs free energy 



G(T ,P ) = E(t)-T S(t) + P V{t); see Ref. 
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for more details. The microstate probability 



Pi(t) in Eq. (!3~Pj) clearly shows the effect of irreversibility and is very different from its 
equilibrium analog p i)eq 

p iieq = exp[/3 (G(T , P ) - E, - P VA]. 

We now provide another demonstration that the two terms in Eq. ( 1231) are identical to 
the generalized heat and work in the following theorem. 

Theorem 5 The heat and work in Eqs. [TU) and ( T771) , and (23) are the same. 
Proof. We first note the identity 

Zi(Vi ~ mEi)d Pi = 0, 

which is nothing but the identity dQ(t) = T(t)dS(t). Using this identity, we obtain another 
identity that follows from Eq. ( j30l ) 

P^iVidpm + A(t)J2^dPi(t) = 0, (32) 

where being related to the Lagrange multipliers, P(t), A(t), etc. must be treated as param- 
eters. We now rewrite dV(t) = £^(t)c^(t) + ErfW(t) and d£(t) = J2Mt)dpi(t) + 
^2iPi(t)dii(t) in dW{t) = P(t)dV(t)+A(t)d£(t) and use the above identity in Eq. fl32} to es- 
tablish that the work in Eqs. (fTTj) reduces to the isentropic form dW(t) = P{t)Y^iPi{t)dVi + 
A(t)J2iPi(t)d£,i(t) . As P(t) and A(t) are parameters, the statistical nature of dW{t) appears 
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in dV(t) and d£(t). As the work expressions are identical, the demonstration also proves 
that the thermodynamic parameters in Eq. (l3Tj) and the statistical fields in Eq. f )24|) are 
the same. ■ 

This proves the two formulations of work to be identical. Apart from providing a check 
for the internal consistency in our approach, this also proves that the first term in Eq. (1231) 
is identical to the heat in Eq. (fit)]) . At equilibrium, A(t) — > A = and the second term 
in Eq. (132]) vanishes. As P(t) — > P at equilibrium, we conclude that d e W(t) = P J2iPidVi, 
see Theorem HI which makes d e W(t) purely isentropic, as expected. 

Let us compute the generalized pressure work W(t) done by the system over some time 
interval t at fixed Pi(t) as V(t) [V(t) = dV(t)/dt] changes from its initial value. It is given 
by the integral 

Wy(t) = EiPiW / Pi (t)Pi(t)V(t)dt = [ P{t)V(t)dt (33) 
Jo Jo 

and differs from the negative of the work W(t) = PoAV(t) (obtained by replacing all the 
quantities in Eq. (|33|) by their medium analogs) done by the medium on the system by the 
dissipative contribution. This result should contrasted with the result obtained by others 



42j . They become identical if and only if the average pressure P(t) is identical with P . 



C. Medium 

The above discussion can be easily extended to the medium (the suffix 5 denotes its 
microstates) with the following results 

dW(t) = -dE _= -52 5 p 5 ^dV = PodV, 
s dV 

dQit) = dE\ = EsEadpa = dQ = -d e Q, 

where all the quantities refer to the medium, except d e Q, and have their standard meaning. 
The analog of Eq. ([29} is 

dQ/dS = (dE/dS^j ~~ = T , (34) 
as expected. We clearly see that 

dW + dW^O 

such as when mechanical equilibrium is not present. In this case, we also have 

dQ + dQ^ 0. 
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D. Irreversible Work and Heat 



We can now identify d{W{t) and diQ{t) : 

d { W(t) = -(dE\ s + dE 
d- x Q(t) = (dE\ VA +dE 

with diQ(t) = d{W(t) from Theorem [TJ Because of the equality, we only needs to compute 
one of them, which we take to be d{W(t), merely because it only involves adiabatic quantities. 
However, the computation of irreversible work requires considering both the system and the 
(working part of the) medium, which makes their computation quite unfeasible in many 
situations because of the very large size of the medium, unless the equations of state of the 
medium are known. On the other hand, the determination of dW(t) and dQ(t) is a much 
easier task computationally as we only deal with the system, a point made several times in 
the literature; see, for example, Jarzynski 

E. External and Internal Variations of dpi(t) 

Let us introduce the natural partition as in Eq. (j2J) for dpf. 

dpiit) = d e pi(t) + diPi(t); 

d e pi(t) is the change due to exchanges with the medium and d\Pi{t) the change due to internal 
dissipation d\Q. Then we have 

d e Q{t) = J2iEid e pi(t), diQ(t) ee Y.iEid-Viit). (35) 

We similarly have 

d e S(t) = J2iVidePi{t), diS(t) =J2iVidiPi{t)- (36) 

We can recast Eq. (JIB]) 

J2i( E i ~ T oVi) d iPi = (T(t) - T )52imdpi 

that acts as a constraint on possible variations d\Pi. The relation d e Q(t) = Tod e S(t) can be 
expressed in terms of d e pi(t) 

~ (3oEi)d e pi = 0. 

We now prove the following theorem. 
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Theorem 6 The identity in Eq. fiEj) is a consequence of the vanishing of average change in 
microstate uncertainty 

which is nothing but the conservation of probability. 

Proof. We express Eq. flS]) in terms of microstate probabilities and use Eq. (1321) to obtain 

E i E l d iPl (t) = P(*)EiPi(*)<ftS + A(t)J2 iPt (t)dUt) - P dV. 
We eliminate the last term using the traditional formulation of the first law to obtain 

$>i(*) [dEi(t) + P(t)dVi + A(t)d£i(t)] = 0. 
In view of Eq. ( 130|) . this is nothing but Y2iPi d Vi = 0- This proves the theorem. ■ 

F. The Adiabatic Theorem 



We now have a clear statement of the generalization of the adiabatic theorem 34] for 
nonequilibrium processes. An adiabatic nonequilibrium process is an isentropic process. 
Such a process also includes the stationary limit, i.e. the steady state of a non-equilibrium 
process. However, the extension goes beyond the conventional notion of an adiabatic process 
commonly dealt with in equilibrium statistical mechanics according to which an adiabatic 
process is one for which d{S = 0, and represents a reversible process in a thermally isolated 
system so that d e Q(t) = 0. One can also have dS = in an irreversible process during which 

diS(t) = -d e S(t) > 0; (37) 

as usual, Eq. ([8]) always remains satisfied. If the system remains in internal equilibrium, 
which may not hold for a stationary state, then we must also have 

diQ{t) = -d c Q{t) = T diS{t) > 0. (38) 

In the following, we will not consider a stationary nonequilibrium state, a case that would 
be taken up elsewhere. 



Theorem 7 In an adiabatic process, the sets of microstates and of their probabilities Pi do 
not change, but d e pi = —diPi ^ for all i. 
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Proof. In terms d e pi and diPi, Eqs. (|37|) and (1381) become 

^ f^iPi - - ^ 7?i<4Pi, (39a) 

i i 

-E^iPi ^4pi- (39b) 

Recognizing that there is only work in the defining relation in Eq.f l23p . which requires Pi not 
to change, we conclude that 

dpi = for Wi 

in an adiabatic process. As diS(t) does not vanish in an irreversible process, d\pi(t) cannot 
vanish. Accordingly, d e pi = —d- x pi ^ for an irreversible adiabatic process. As p^s do not 
change, no microstate can appear or disappear. This proves the theorem. ■ 



G. Ideal Gas: An Illustrative Example 

Consider, as an example, an ideal gas in a cuboid with its axis along the x-axis. Its 
length a(t) along the x-axis is controlled by a frictionless piston at one end. All other walls 
are assumed rigid. Let the dimensions b and c, respectively, along the y- and z- axes remain 
constant. As the particles are ideal, we can focus on single particle energy levels given by 

' 2m a 2 (t) b 2 c 2 
with the corresponding wavef unction given by 

0n(r(t)) = \ -7WT- sin(-^rx) sin(^|y) sin(— z), 
]l a(t)bc Tra[t) -Kb ttc 

which form the normalized complete basis set. Let p n (t) denote the probability of having 
gas particles in a given eigenstate n. Then the energy per particle of the gas is given by 

E(t) = J2Pn(t)E n (t), 
n 

and the pressure on the moving piston is given by 

„. . 1 . dE . h 2 ir 2 . . n? 2 trl , . _ , 

p « = s(-sw) = imt) 5>«>4) = W) m) ~ E ° l 

where 

/ x h 2 n 2 ,nl 

E » = EM*) 1 ^(jf + ^). 

n 
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The probabilities and the energy eigenvalues change in time as a function of a(t). The 
change in the energy E(t) comes from changes in p n (t) and from the variations in a(t). The 
two contributions determine heat and work, respectively. 

If the system of ideal gas is isolated, its energy, volume and the number of particles remain 
constant. If the gas is originally not in equilibrium, it will eventually reach equilibrium in 
which its entropy must increase. This requires the introduction of some internal variables 
even in this system whose variation will give rise to entropy generation by causing internal 
variations dip n (t) in p n (t). Here, we will assume a single internal variable Its real 

nature is not relevant for our discussion. What is relevant is that the variation in £(£) is 
accompanied by changes d(p n (t) occurring within the isolated system. According to our 
identification of heat with changes in p n (t), these variations must be associated with heat, 
which in this case will be associated with irreversible heat diQ{t). 

1. Isothermal Expansion 

Let us first consider an isothermal expansion of this gas so that the temperature of the 
gas remains constant and equal to that of the medium T . During expansion, energy is 
pumped into the gas isothermally from outside so that not only E(t) remains constant, but 
also keeps P(t)V(t). The pumping of energy will result in the change d e p n (t). This will 
determine d e Q(t) = T d e S(t). In addition, particles may undergo transitions among various 
energy levels, as discussed above, without any external energy input, which will determine 
the change dip n (t). This will determine dsQ(t) = T diS(t), and consequently d\W(i). Thus, 

[Pit) - P ]dV(t) + A(t)d£(t) = T dS(t) - d e Q(t), 

which allows us to determine the irreversible work in terms of measurable quantities. Such 
a calculation will not be possible in the traditional formulation of the first law. 

2. Adiabatic Expansion 

In a nonequilibrium adiabatic process, we have d\W(t) = — d c Q(t) so the heat exchange 
\d e Q(t)\ = To \d e S(t) \ is converted into the irreversible work in this process. We can use this 
to determine the work d[W^(t) due to the internal variable 

A(t)d£(t) = -d e Q(t) - (P(t) - P )dV > 0. 
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The identification d)W(t) = — d c Q(t) and the calculation of A(t)d£(t) cannot be done in 
the traditional formulation of the first law. 



VI. CLAUSIUS EQUALITY 

It follows from Eq. ( Tit)]) that dQ(t) /T(t) is nothing but the exact differential dS(t) so 
that 

j) dQ(t) /T(t) = for any cyclic process. (40) 

We will call it the Clausius equality. The equality should not be interpreted as the absence of 
irreversibility; see Eq. @. It is only because of the use of the generalized heat dQ(t) in place 
of d e Q(t that the Clausius inequality has become an equality. Using diS(t) = dS(t) —d e S(t), 
we obtain the original Clausius inequality for a cyclic process taking time r: 

N(t, t) = fdiS{t) = -§d c Q{t)/T > 0, (41) 

which is the second law for a cyclic process and represents the irreversible entropy generated 
in a cycle; see Eq. (j3J). The quantity N(t) is the uncompensated transformation of Clausius 



6| that, as we have just discovered, is directly related to diS(t) [50J]; in contrast, N (t,r 



N Q {t,r) = §d- 1 Q{t)/T{t) (42) 

is determined by the uncompensated heat d\Q(t) and represents a different quantity. In terms 
of the two heats, we have 

§d i Q{t)/T{t) = -§d e Q{t)/T{t) > 0, (43) 

which results in a new Clausius inequality 

§d e Q(t)/T{t) < 0; (44) 

compare with the original Clausius inequality in Eq. (j2J). The contribution No(t) in Eq. 
f|4*3|) can be thought of as a new uncompensated transformation, different from N(t) [6j|. Our 
formulation has allowed us to identify the two transformations 

N(t,r) ee faS® = §<kQ{t)lT Q -§T{t)dS(t)lT Q , N (t,r) = §<kQ{t)/T{t) = §W(t) /T(t), 

which provides a way of computing them using nonequilibrium equations of state of £ and 
E. 
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Other quantities also appear as equalities. Let us consider work, which appears as 
dW(t) = — \dE{t) — T(t)dS{t)}. For a not too small a change, we have 



AW(t) 



AF(t) = A [E(t) - T(t)S(t)] . 



(45) 



The quantity F(t) in the brackets should not be confused with the Helmholtz free energy 
F(t) = E(t) — T S(t). The above equalities should be contrasted with the traditional 
inequalities d e W{t) < -dF{t) or A c W{t) < —AF(t). 

VII. APPLICATIONS 

We now consider two simple applications of the new approach. 

A. Free Expansion 

Consider the example of free expansion in a set up in which gas in one chamber is 
separated from another empty chamber of identical volume V° by a partition; the latter 
is held stationary by some mechanism. Both chambers form an isolated system so that 
not only d e Q = but also that the expansion occurs at constant energy E°. Therefore, 
d e W = 0. Therefore, the traditional formulation of the first law is of no use in obtaining any 
useful information about the nature of irreversibility. The initial pressure and temperature 
of the gas are P° and T°. As soon as the mechanism is removed, the partition becomes 
free to move and the gas expands into the other chamber against a zero pressure. The 
pressure difference AP(t) = P(t), where P(t) is the instantaneous pressure of the expanding 
gas occupying a volume V(t); the latter can be recorded as a function of time. Using 
the recorded V(t), the temperature T{t) and P(t) of the gas are then obtained from the 
instantaneous equation of state of the gas and the condition E (t) = E°. The resulting 
irreversible work is djW{t) = P{t)dV{t) = d^Q(t) = T(t)diS(t), where we have used Eq. 



where the suffix E is used to represent that E is constant in the process, and where S\ E 
represents the initial entropy. Neither the initial state nor the state at t is required to be an 



f TT6b]) with d e S(t) = 0. We thus find 




o 

E ■ 
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equilibrium state in the above calculation, which follows the recorded non-equilibrium path 
given by V(t) at constant E. 

Assuming the initial and final states to be equilibrium states, we have 

A { S{t oo) = S cq (E°, 2V°, N) - S cq (E°, V°, N). 

For an ideal gas, we evidently obtain A[S(t — > oo) = iVln2 between the two equilibrium 
states, as expected. 

The work done by and the heat generated within the system are given by 

r v(t) 

AW(t) = AiQ(t) = AiW(t) = / P(t)dV(t), 

Jv° 

and can be calculated for the expansion profile given by V(t). 
B. Relative Motion and the Resulting Friction 



Let us consider a different situation, in which 
the medium with some velocity V(t) given by 



;he system is also moving with respect to 
M 



0(t)V(t) = - (dS(t)/dP(t)) E{t)y{t) 

in which P(t) and E(t) are the momentum and internal energy of the system in the lab 
frame (in which So, and to a very good approximation £ are at rest). Such a motion arises 
during sudden mixing of fluids or in a Couette flow, and eventually stops due to deceleration 
caused by generated frictional forces as a result of the relative motion; we do not consider 
any external force causing this motion. The irreversible work done by the system against 
friction is 

dfW(t) = -V(t)-F fr dt > 0, 

where Ff r is the resulting frictional force opposing the motion and results in dP(t) = Ff T dt. 
It is clear that dfW(t) vanishes as the motion ceases. This dissipative work is in addition 
to the irreversible work caused by any pressure difference and any affinity as before: 



diW(t) = (P(t) - P )dV{t) + A(t)d£(t) - V(t)-dP(t). 
From Eq. (jT6l) and Eq. ([8]), we then find 

T(t)ckS(t) = [(P(t) - P )dV(t) + A(t)d£(t) - V{t)(t)-dP(t) - (T(t) - T )d e S{t). 
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(46) 



The same conclusion can also be obtained by considering the Gibbs fundamental relation 
for this case. Using the above definition of the velocity, it is easy to see that 

dE(t) = T(t)dS(t) + V(t)-dP(t) - P{t)dV{t) - A{t)d£(t)- (47) 

see Ref. [8j. The last three terms, being isentropic contribution to dE(t) represent dW(t), 
while the first term represents dQ(t). We now follow Ref. [8] and rewrite 

dW(t) = P dV(t) - Vo-dP(t) + (P(t) - P )dV(t) - (V(t) - Vo)-dP(t), (48) 

with the first two terms representing d e W(t), and the last two terms d[W(t). Here, Vo 
represents the equilibrium pressure velocity. However, the equilibrium value of V(t) is 
V = 0, which then yields djWit) given in Eq. f|46|) . The exchange work d e W(t) = P dV(t) 
remains unchanged and says nothing about the presence of internal friction terms that 
appears in our approach. 

VIII. INCLUSION OF OTHER STATE VARIABLES 

Let us now extend the discussion to include other extensive quantities such as the flow 
of matter, the electric interactions, chemical reactions, etc. For specificity, we focus on 
chemical reactions, which we assume to be described by a single extent of reaction £(t). The 
corresponding affinity for the system is given by A(t), while that for the medium is given by 
Aq = 0. We assume another observable X such as the number of solvent in a binary mixture. 
The corresponding chemical potential is /i(t) for the system and hq for the medium. The 
work is now 

dW(t) = P(t)dV(t) - n(t)dX + A(t)d£(t). (49) 
The Gibbs fundamental relation for £ is given by 

dE(t) = T(t)dS(t) - P{t)dV{t) + n{t)dX - A(t)d£(t), (50) 

while the first law for it takes the form dE(t) = dQ(t) — dW(t). Rewriting dQ(t) = dE(t) + 
P(t)dV(t) - n{t)dX + A(t)d£(t) as dQ(t) = dE(t) + P Q dV(t) - ^ dX + (P(t) - P )dV(t) - 
(n(t) — fio)dX + A(t)d^(t), we can identify [8] d e Q(t) with the first three terms 

d c Q(t) = dE(t) + P Q dV(t) - HodX. 
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The last two terms represent d e W(t): 

d e W(t) = P dV(t) - HodX, 

which shows that the second term will appear in the traditional formulation of the first law, 
but not the associated thermodynamic force fiQ — fi(t); see below. Similarly, d\Q(i){t) = 
diW{t) represents the dissipation given by 

diW(t) = (P(t) - P )dV(t) - (ji(t) - fi )dX + A(t)d£(t); (51) 

see Eq. (JSJ). According to the second law, each contribution in d\W(i) is non-negative. The 
irreversible entropy generation is given by 

TodiSit) = [T - T(t)]dS(t) + [P(t) - P )dV(t) + A(t)d£(t) + (/i - fi{t))dX(t). (52) 

The extension to arbitrary number of observables and internal variables is trivial js]. 



IX. A CLOSED SYSTEM PERFORMING WORK ON AN EXTERNAL OBJECT 

We now consider our isolated system to consist of the previous medium E(T , Pq, A ), the 
system E[T(t), P(t), A(t), fi(t)} and an external object S cxt [T , Pq, Aq, /i ext (t)]; the latter is 
completely disconnected from the medium but can only exchange the type of work with the 
system involving the extensive variable X(t) considered in the previous section. However, 
there are some important differences. First, we do not limit the external object to be 
extremely large compared to the system. Accordingly, its chemical potential /z e xt(£) can 
change in time. Just for convenience, we assign to it To, Pq, Aq of the medium, even though 
it is isolated from the latter and cannot exchange heat and volume and internal variable 
work with the system and the medium. It is a classic problem that has been extensively 
studied in which the external object is thermally insulated form the system and the medium; 
see for example Ref. [3J, see Sec. 20 there]. It is also studied by the stochastic trajectory 
approach; see for example, Refs. 0] • Here, we will closely follow Landau and Lifshitz 
3J, see Sec. 20 there]. As the external object only performs work, its entropy S ext must not 
change. We will denote the totality of E and E ext by E and all its quantities by an additional 
bar. For the isolated system E, we have AS (t) = A i S , (t) = A S(t) + S cxt (t) + S(t) = 
S(t) + S(t) > during a process involving (not necessarily small) change within the 
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isolated system E . As usual, all exchange quantities between various pairs of the three 
parts cancel out for So; only irreversible quantities survive such as A[So{t), which in our 
approach is associated with the entropy generation within S. Accordingly, we will use the 
notation A i S'(t) for A i S'o(t) below and write 

AiS(t) = A;S (Q) (t) + A ; 5 (v) (t) + AiS^it) + A^ (X \t) > 0. (53) 

The first three terms on the right in the first equation denote the cumulative irreversible 
entropy changes due to exchanges between the system and the medium and the last term is 
the cumulative irreversible entropy change due to the exchange between the system and the 
external object. The work done A e R(t) by the external object on the system is 

A c R(t) = A [E(t) - T S(t) + P V(t)} + T Q AiS(t). (54) 

In terms of A c W(t) = —A c R(t), we have the following inequality 

A e W(t) < AG(T , P , t)=A [E(t) - T S(t) + P V(t)] (55) 

where AG(Tq, Po,t) is the change in the Gibbs free energy of the system. It should be 
stressed that the above inequality does not assume the establishment of internal equilibrium 
in the system. If we are considering spontaneous changes occurring within the system in the 
medium S(T , P , A ), we must set the external work d e R = above by setting dX(t) = 0. 
We then find 

dG(T ,P ,t) = -T d i S(t)<0 (56) 

during spontaneous relaxation. Thus, assigning the irreversibility to the system leads to the 
expected behavior of the time-dependent Gibbs free energy of the system. 

The corresponding work done by the system with respect to the variable X(t) is dWx{t) = 
—fi(t)dX(t), while the work done by the external object is d e R(t) = /i ext dX(t), with the 
result diW x (t) = [/i ext (t) - //(*)] > 0. We then find 

TodiSit) = [T -T(t)}dS(t) + [P(t)-P }dV(t)+A(t)d^t) + [fi ext (t) -fi(t)]dX(t) > 0, (57) 

in which the first three terms refer to irreversible entropy generation within the system and 
the last term refers to the generation within E. Each term on the right must be non- negative; 
compare with Eq. ( TT5|) . Comparison with Eq. (15"TT) shows that the last three terms above 
represent diW(t), except that /z is replaced by [i ext . The additional thermodynamic force 
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A*ext(^) — >u(t) vanishes when there is equilibrium established between the system and the 
external object; let the equilibrium value be given by ~p . The new force generates the 
irreversible entropy diS (t) due to X. Following the arguments given in Ref. [8[, see Sec. 
XI there, for the case of thermal equilibration and which can be adapted here verbatim, we 
find that the irreversible entropies for £ and E ex t are given, respectively, by 

TodiS^it) = [7J -^(t)]rfX(t),To^2(t) = [p -^(t)]dX ext (t), (58) 

with dX(t) + dX eXk (t) = 0. As the entropy of the external object remains constant, there 
must be a corresponding change d e S^(t) = —diS^(t) as it performs work, which requires 
d e Pk,ext(t) = —diPk,cxt(t) for the k-th microstate of the object. 

If the object is extremely large compared to the system so that the coupling with the 
system will not appreciably affect fi ext , the latter can be taken as a constant equal to fio = /Z - 
In this situation, the entire contribution d^S (t) must be treated as diS^(t) assigned to 
the system, as there is no force in the object to generate irreversibility (diS ext (t) = 0). 

The difference between our approach and the stochastic trajectory approach should be 
pointed out. In the latter approach, T = T Q and P = Pq for the system. Also, as no internal 
variable is considered, we will set A = A = above. Thus, the only irreversibility is due 
to the possibility of work by the external object. The entropy S(t) is given by the Gibbs 
formulation; see Eq. (I2ip . and its changes dS(t) includes the irreversible entropy changes 
diS(t) = diSW(t). This is not the case in the stochastic trajectory approach in which the 
irreversible entropy AiS^it) is assigned to the medium [42|]. 



X. SUMMARY AND DISCUSSION 



A. Summary 

The traditional formulation of the first law has two basic deficiencies. The first one is that 
even though the first law is about the change in the energy of the system, it contains the 
fields of the medium and not of the system. The second, and more dramatic conceptually, 
deficiency is that it remains satisfied even if the process violated the second law. This is why 
both laws are needed in nonequilibrium thermodynamics. Despite the concern expressed by 
Sect. 5.12], we have shown that the definition of the generalized heat dQ(t), 



Kestin 



12 



which includes its irreversible component d]Q(t), follows uniquely from the unique choice of 
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dW(t) resulting from the second law, see Eq. (|T9|) . and also from the statistical approach, 
see Eq. fl25|) . It general, we have the identity d\W(i) = d[Q(t) between their irreversible 
components; see Eq. ([H]). These results do not require the system to be in internal equilib- 
rium. The situation becomes very illuminating when the system is in internal equilibrium. 
In this case, the generalized heat formulation gives rise to not only the simple relation- 
ship dQ{t) = T(t)dS(t), but also the extension of the adiabatic theorem for nonequilibrium 
states; see Theorem [7J Another remarkable consequence is that in terms of generalized 
dQ(t), the Clausius equality ( |40|) is always maintained, in contrast to the inequality ( H3|) 
in the traditional approach. Our generalized formulation brings about a symmetry between 
the system and its surrounding medium, see Eq. ( II 9p . which is absent in the traditional 
approach using external fields. The first law becomes identical to the second law so that 
we only deal with equalities and a single law; the equalities are easier to deal with than the 
inequalities that result in the traditional approach; see for example, Eq. (H5|) or Eq. (J7J). 

The introduction of microstate probabilities and the partition of dpi(t) into d e pi(t) and 
diPi(t) provide a very convenient and useful statistical representation of exchange and ir- 
reversible forms of heat and entropy and of other quantities, and which allows us to draw 
some useful conclusions. For example, it follows directly from Theorem ([7]) that pure work 
cannot change the entropy. This was used to argue that S ex t must be a constant in Sec. IIXI 
We have found that the microstate probability Pi(t) is very different from its equilibrium 
analog Pi, e q- Several examples have been given to show the usefulness of the new definition 
of heat and work. Below we list some of the benefits of the new formulation. 

B. Benefits of the Generalized Formulation 

1. It is a new way to write the first law, which shows that the first and the second laws 
are no longer independent, a property now shared by all irreversible processes and not 
just by reversible processes. This conclusion must be contrasted with the well-known 
result found in all text books and the literature that the two laws are independent for 
irreversible processes. 

2. Work and heat can be uniquely distinguished so that each has a clear physical signifi- 
cance, but no such clear identification is possible in the traditional formulation of the 
first law; see Theorems ([3]) and @. 
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3. Work and heat no longer are restricted to flow through the boundary. Thus, the new 
formulation can easily accommodate interactions described by fields (electromagnetic, 



gravitational, etc.) [12 1. 



4. The entire thermodynamics and stability considerations for any body can be expressed 
in terms of the variables associated with the body alone at each instant. We can use the 
instantaneous equations of state of the body alone for thermodynamic computation. 
Indeed, it may be the only way to study a non-equilibrium system in some cases as 
is clearly shown by the example of free expansion in Sect. I VIII Here, the traditional 
formulation of the first law is of no use to determine d{W{t) and diS(t) as is well 
known. 

5. We deal with equalities and not inequalities. In addition, we only deal with dS(t), 
which as a dependent variable is expressed in terms of state variables Z(t); there is no 
need to consider its parts d e S(t) and diS(t) separately, which depend on the medium; 
see (3) in Sect. Ill PL The generalized heat dQ(t) and work dW(t) differ from d c Q(t) 
and d e W(t), respectively, by the same contribution diQ{t) = d\W{t) = diWoit). 

6. The use of the Gibbs fundamental relation makes it almost trivial to identify the 
irreversible contributions. Thus, the new formulation does not have the limitation 
inherent in the traditional formulation; see (7) in Sect. Ill Dl The determination of 
diWo(t) is straightforward by measuring P(t) and P . Adding this to d e Q(t) then 
allows us to determine dQ(t). 

7. The most useful aspect of the generalized approach is that the work dW(t) does not 
depend on the amount and nature of dissipation going on during a process. It is 
determined by the instantaneous pressure of the system, which itself is a derivative of 
its instantaneous entropy, a state variable. We have already seen this earlier in Eq. 
fll2p . Whether friction between the piston and the walls is present or not, d\W(i) is 
always given by AP(t)dV(t) so that dW{t) = P{t)dV{t). 

8. Once dQ(t) is known, dS(t) = dQ(t)/T(t) is also known or vice versa as if we are 
dealing with a system in equilibrium, even though we are dealing with an out-of- 
equilibrium system in internal equilibrium. In the latter situation, the analysis is 
considerably simplified as we are always dealing with the instantaneous properties such 
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as the equation of state of the system even though the system is out of equilibrium. 
We do not need to, or may not, know these properties for E. 

The new formulation has many other desirable properties. The generalized heat dQ(t) and 
work dW(t) differ from d e Q(t) and d e W(t), respectively, by the same contribution diQ{t) = 
djW(t). The determination of d{W(t) is straight forward in terms of thermodynamic forces 
P(t) — Po,A(t), etc. Adding this to d e Q(t) then allows us to determine dQ(t) from which 
dS(t) can be determined. 
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